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Abstract 

We study the damping of low-lying collective excitations of condensates in a weakly 
interacting Bose gas model within the framework of imaginary time path integral. A 
general expression of the damping rate has been obtained in the low momentum limit 
for both the very low temperature regime and the higher temperature regime. For the 
latter, the result is new and applicable to recent experiments. Theoretical predictions 
for the damping rate are compared with the experimental values. 
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The recent realization of Bose- Einstein condensation (BEC) in dilute atomic vapors |1|, 
has opened the door to experimentally study weakly interacting dilute quantum gases for 
which microscopic theories have been well studied for decades. In particular, it has been 
shown that the frequencies of the low-lying collective excitations of the condensates P, Q 
agree excellently with theoretical predictions based on a mean-field theory . The damping 
in such collective modes has been experimentally discovered by Jin et al and Mewes et al 
A very recent experiment has extended the study of low-lying collective excitations 
of condensates to include higher temperatures in which the damping of these excitations 
exhibits dramatic temperature dependence. So far there is no theoretical prediction for it in 
such a temperature regime. 

This paper aims to provide a simple way of calculating the damping rate of collective 
excitations in the above temperature regime. The experimental setup we are considering is a 
dilute gas of atoms confined in a trap potential. Recently, Chou, Yang and Yu []7| argued that 
the local density approximation can be applied into the present problem. In other words, 
the trap potential may be treated as a slowly varying external potential and enters into 
the theory only as a modification to the chemical potential. To simplify the theory, let us 
assume that the trap is absent in deriving the damping rate since we are only interesting in 
finding the leading term for it. In the following, we shall study a model of weakly interacting 
bosonic particles in an infinite free space using the imaginary time path integral method 
which allows us to take into account excitation collision processes that are believed to be the 
causal of the damping. 

In the low momentum and low temperature limit, one usually finds an effective field 
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theory in which only slowly varying degrees of freedom appear explicitly with interactions 
that include the effects of fast varying fields that have been integrated out first One 
version of such an effective field theory has been established for a non-ideal Bose gas in the 
low temperature regime by Popov [P, [T^. The theory treats the system effectively as that 



of quasiparticles described by two slowly varying real scalar fields and cr(x). Here, the 
cr(x) field describes the density fluctuations, written as a{x) = n{x) — no with uq the density 
of the ground state. The field is the Goldstone field in this theory appearing originally 
in the phase of original particle fields. Notice that we have written the four- dimensional 
Euclidean spacetime as a; = (x, r) with r = it denoting the Euclidean "time" . According to 
Popov, the Euclidean action for such a system can be written hj {h = 1 henceforth) 

+-^a^( ) cr(x)(V0(x)- V0(x)) | 

2(9n2^ 8mno 2m J ^ ^ 

where m is the particle mass and p is the pressure of the system given as the function of 
the chemical potential /i and the particle density n. Notice that and cr(x) are periodic 
in "time" r with period f3 = 1/(A;bT). If v = (l/m)V0 is identified with the phonon 
velocity field, one can easily verify that the action (|l]) corresponds to a Hamiltonian that is 



one form of Landau-Khalatnikov hydrodynamic Hamiltonian |Tl|. In the low temperature 
regime such that the non-condensate fraction of particles is very small, the expression of the 
pressure p{fi, n) can be well approximated by that of T = 0. For a weakly interacting dilute 
gas, it is given by [|T2| 

^0 2 

p = (in — —n 

where to = ^ characterizes the interaction of particles with the s-wave scattering length a. 
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It follows that 
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For simplicity, we write and a fields into a real scalar doublet as $ = (0, a) . The 
Green's functions (or propagators) are defined by the matrix 



G{x-x') = {T{^{x)^{x')}) 



I[Ux d(f){x)da{x)]<^>{x)<!>{x')e^^'^'''^ 
m.d<P{x)da{x)]e'^lM 



with T denoting a time-ordered product. The Fourier transform of G{x) is defined through 



(3) 



where the notation k = {k,iui,) is understood and Matsubara frequencies tOy = 27rz///3 
(z/ = 0, ±1, ±2, ■ ■ ■). The quadratic part of the action 5* can be written into the form of 
'S'quad = / d'^xd'^x'^{x)^{x,x')^{x') and the free Green's function G'o(^) is equal to the 
inverse of the matrix Thus, we find 
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It follows that 



I to + A;V(4mno) -u^ \ 
tUjy {no/m)k'^ 



(5) 



where the spectrum 



(6) 



with c = ytQUo/m. The propagators Go are represented by Feynman diagrams in Fig. |I[ 
Further, the cubic term of the action (|I|) describes the interaction of three excitations that 



are known as phonons in the low momentum region, giving rise to a vertex of 

(ki ■ k,) 



5^(ki + k2 + k3)(5, 
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represented by the last diagram of Fig. |I[ Obviously, it can be treated as perturbation in 
the low momentum limit. 

The spectrum of collective excitations is given by the poles of the exact Green's function 
[1T3| G{k) = Go{k) + Go{k)Il{k)G{k) where n(fc) denotes the matrix of the self-energy parts. 
It follows that the spectrum is determined by 

det G-\k) = det{Go\k) - U{k)) = 0. (7) 

Now that we are working in the imaginary time formalism, let us make the analytical contin- 
uation ^ u + ir] {r] = 0^) after the Matsubara frequency sum and write u = Reu — i'y(k) 



with 7 denoting the damping rate. Then, keeping only terms up to one-loop order we 
find from Eqs. @ and (|^ that 



7(k) ^ 



2Reu 



(to + -; )lmn<^0(k, u + ir])-\ Imn^^(k, u + irj) 



Rell^^(k,u + irj) 



with the understanding that all cj's are replaced with e(k) after the analytical continuation. 
Up to one-loop order, there are six diagrams (see Fig. ^ for the self-energy matrix n(k, uj) 
but only the last five (b)-(f) contribute to the damping. After collecting contributions from 
all related diagrams, we have 

7(k)=7i(k)+72(k), (9) 

with 

7i(k) = ^/rf3A:'5(6(k)-6(k')-6(k-k'))[/(6(k'))-/(-6(k-k'))]x 



(k - k'Yik ■ k')'e(k')e(k) (k • k')(k • (k - k'))e(k) 



2mnoFfc'^e(k - k') 
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and 



k\k' ■ (k - k'))'e(k')e(k - k') (k' • (k - k'))(k ■ k')e(k') 1 
4mno/i;'^(k - k')2e(k) mnok''^ J 



72(k) = ^/rf3F<5(e(k) + e(k')-e(k + k'))[/(e(k'))-/(e(k + k'))]x 



(10) 
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where /(e) = l/(exp(/3e) — 1). 

If T = 0, one can easily verify that for small k {ck ^ noto) 

3k^ 



(k + k')2 



(11) 



7T=o(k) 



(12) 



which is the well-known Beliaev's result [15|. The product rao^o characterizes the strength of 



particle interactions. 

For T 7^ and small k such that ck <^ ksT and ck -C noto, we find that the damping 
rate to the lowest order in k is determined by 



7(k) 



ckk^ ^^ Jiptp ^ 
IQ-nmnQkBT ksT' 



(13) 



where k^ = ^JmriQtQ and the function ^{x) is given by 

2e 3 



xix-^ Jo 
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If the limit ^ is taken, Eq. (|TB[) reduces to the familiar form 



7(k) 



37T^k{kBTY 

40mnoC^ 
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(14) 



which was given by Hohenberg and Martin |]T6|. Though, we shall see in the following 



paragraph that Eq. ([1^) is invalid in the recent experimental temperature regime. In fact, 
both expressions ([T^ ) and were already derived by Popov in Ref. from the effective 
action (|l|). But, to the best of my knowledge, the general expression (|13|) is first obtained 
here. It is valid as long as the system is in the low temperature region such that the number 
of particles in the excited states is much less than that in the condensate. 

The damping of collective excitations in BEC has been measured in the dilute atomic 
vapors of both ^''Rb and sodium at JILA and MIT, respectively [Q, ^, |]. For the experiment 
of ^^Rb at JILA, the condensates are produced in a trap with frequencies of z/,, = 129Hz 
radially and = 365Hz axially. For a typical condensate of Nbec = 4500 atoms, the 
density of a condensate can be estimated as no ~ 1.4 x lO^^cm^^ Further, the scattering 
length for a ^''Rb vapor may be taken roughly as a = lOSBohr |T3]. It turns out that the 



typical interaction energy noto ~ 57nK. Therefore, the damping rate is expected to be 
given by Eq. ([TsD instead of Eq. ([l^) in the temperature region (30, 300)nK where it was 
measured (see Figs. 1 and 3 of Ref. @). We also checked that for the lowest excitation mode 
(m = 2) both ck/lksT) and ck/lrioto) are less than 0.2 for any T above 50nk. Fig. |^ shows 
theoretical predictions for the damping rate of collective excitations based on Eq. (|13D in 
comparison with the experimental data from Ref. [^. The theoretical curve (6) seems to fit 
the experiment well apart from that the theoretical values are varying with the density of 
condensates approximately according to 7 oc However, a Bose condensate in a trap 

is essentially inhomogeneous and we will not try to determine a proper average density uq 
here. Next, let us check whether the above experiment is in the low temperature region. 
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For T < lOOnk, we read off from Fig. 1 of Ref. |^ that the corresponding T' < 0.6 (where 
T' = T/Tc) and the condensate fraction Nbec/N is greater than or about 80%. Hence, the 
higher order corrections to Eq. (|T3|) due to finite temperature effects only can be estimated 
from the calculation of Feynman diagrams less than (1 - Nbec/N^ = 4%. On the other 
hand side, Ref. ^ shows that more than 50% of particles reside in noncondensates when 
T' > 0.8 and correspondingly T > 175nK roughly. This means that the damping rate shall 
be no longer determined by Eq. ( |T3|) in that temperature regime. I do not know yet whether 
there exists any theoretical prediction for it. 

In Fig. ^, we plot the damping rate of the collective excitation of frequency 30Hz for 
the sodium gas system of Mewes et al 0] using the s-wave scattering length a = 65Bohr 
20| . The decay time of 250(40)ms was found experimentally when a nearly pure condensate 



was formed at the temperature T ^ O.ST^ [pil]. The Bose-Einstein transition temperature 
is determined theoretically according to Tc = |^(A^/ 1.202)^/^ |^ ^ with lj the geometric 
mean of the harmonic trap frequencies Co = (cUx^^yiUzY^'^ . For that the trap used in Ref. 
is of frequencies of 250Hz radially and 19Hz axially and typically holds a total number 
N ^ 5 X 10^ of atoms when the decay time was measured, we have Tc ~ 800nK. For a 
typical condensate density no = 3 x lO^^cm"'^ pO[, we can read off from the curve (c) of 



Fig. ^ that 7 ~ 4.4s~^ for temperature T = 400nK. This damping rate corresponds to a 
decay time of about 230ms that agrees very well with the experimental value. Also, since at 
such a temperature the condensate fraction of atoms is around 90%, higher order corrections 
to 7 due to the finite temperature effect alone are quite small (< 1%). Also, we checked that 
both ck/{kBT) and c/c/(noto) are less than one percent. That is to say, the system is well in 
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the low momentum region. 

In conclusion, this paper analytically calculates the damping rate of collective excita- 
tions for a dilute Bose gas model in a temperature regime where theoretical predictions did 
not exist previously. Although the model has ignored the contribution of the trap and has 
required that the condensate be homogeneous, it produces results for the damping of col- 
lective excitations that are in good agreement with the experiments. Our study also reveals 
that the damping is due to the three-phonon collision processes at finite temperature. This 
paper serves a very preliminary study of this subject and fully understanding it requires a 
complete analysis to take into account the contribution of the trap and the inhomogeneity 
of the system. 
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I am grateful to E. Cornell and W. Ketterle for very helpful instructions on the experiments, 
to M. R. Matthews for sending me the experimental data, to G. Ordonez for helpful con- 
versations, and to Professor E. C. G. Sudarshan for illuminating discussions on BEG. This 
work is supported by NSF grant PHY-9511632 and the Robert A. Welch Foundation. 
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Figure 1: Diagrams for propagators and vertex. 
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(d) (e) (f) 



Figure 2: One-loop diagrams for the self-energy matrix Il{k). 
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T(nK) 



Figure 3: The damping rate of collective excitations in a ^ Rb atomic gas. Three solid lines 
indicate theoretical predictions from Eq. (|TB]). For those, the excitation frequency is taken as 
uj/2tt = lAur = 180.6Hz corresponding to the mode m = 2 and the condensate densities are: 
(a) no = 1.0 ■ lO^^cm-3; (b) uq = 2.0 ■ lO^^cm^^; and (c) uq = 3.0 ■ W^cm'^. Two discrete 
curves are replotted from the data of Ref. 0. (The experimental data is the courtesy of M. 
R. Matthews pf.) 
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Figure 4: The damping rate of collective excitations in a sodium atomic gas. All solid lines 
are plotted according to Eq. (pTSl). The excitation frequency is 30Hz and the condensate 
densities are: (a) uq = 1.0 ■ lO^^cm"^; (6) uq = 2.0 ■ lO^^cm"^; (c) hq = 3.0 ■ lO^^cm^^; and 
(rf) no = 4.0 ■ lO^^cm"^ 
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